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SHELL  THEORY  FROM  THE  STANDPOINT 
OF  FINITE  ELASTICITY 


P.  M.  NAGHDI 

UNIVERSITY  OF  CALIFORNIA,  BERKELEY 


ABSTRACT 

Thir  is  an  account  of  the  nonlinear  theory  of  thin  shells  from  the  stand  - 
point  "f  finite  nr t i c i t y , togcUc-r  with  a brief  summary  of  some  related  recent 
r-  -ircli* s on  the  subject..  The  development  of  the  basic  theory,  as  well  as  n 
liseussUn  of  constitutive  equations  for  elastic  shells,  are  presented  via  a 
lira>ct  approach  on  the  basis  uf  a continuum  model  known  as  a Cos  serai  surface 
rath*,  r than  from  the  three-dimensional  equations  of  nonlinear  elasticity. 

1.  1 INTRODUCTION . GENERAL  BACKGROUND. 


Tii*1  development  of  a complete*  two-dimensional  theory  of  tliin  elastic  shells 
from  tii-  three-dimensional  equations  of  classical  continuum  mechanics  is,  in 
■ •"ivni,  a difficult  problem.  Most,  of  the  difficulties  arise  in  the  development 
< f ‘.h*-  constitutive  equations  and  remain  ev<n  when  the  deformation  is  small . The 
nature  of  these  difficulties,  in  n derivation  from  three-dimensional  equations  has 
b*  n • iai  orated  upon  previously  by  Naghdi  (.1,  Lection  1|.  Because  of  tie-  dil'- 
f i • 1 1 1 1 i • just,  referred  to,  when  dealing  with  motions  and  deformations  of  thin 
sh  li-li’--  bodies,  it  has  be*  n customary  to  employ  a variety  of  approximate 
l-r  • iur  i n t h-  • i « *r  i v • • 1 ion  of  elastic  shell  theories  from  1 he  thro*  -1  imenr  ion  « I 
ii-iti  n.  . 1 r i : • t • d o 1*  adopting  such  n procedure,  w-  approach  the  subj'-ei.  h>  n 

I r*  •**.  ■ • 1 1 < ' t * i r point  of  view,  namely  via  the  theory  of  a 1 • r • * 1.  (or  a d i r ct  ■ d 

. i r* f ■ 1 • ’ • ' . which  is  bas*'d  on  a continuum  model  comprising,  a material  surface  in 
l’  e ' i b an  V-pace  with  n deformable  vector  field  --  called  a director  --  attached 
• • • r..  point  of  the  surface.  It  should  be  emphasized  that  a Cosserat.  surface 

is  not  just  a two-dimensional  surface  but  one  which  is  endowed  with  structure. 

A genera i background  and  motivation  for  use  of  a direct  approach  in  shell 
tii-  ry,  as.  veil  as  a further  description  of  a Cossernt  surface  together  with  som< 
historic  1 background  on  Hie  subject  can  be  found  in  flj.^  Here,  we  limit  our- 
. t..  only  ••  erief  historical  sketch  on  the  basic  theory  of  'osserat  surface 

a.  con c pt  of  'directed*  or  'oriented*  media  originat'd  in  the  work  of  Du  hem  I .J 
ii  i a first,  systematic  development  of  theories  of  oriented  media  in  one,  two  and 
Hire  -dim* ’ns ions  ( th-  first  two  being  motivated  by  rods  and  sheLls)  was  carried 

so,  s-  ferine-  may  be  made  to  Naghdi  [2j  which  contains  a rapid  exposition  of 
th  ry  of  a boss- rat  surface  and  includes  a discussion  of  its  relevance  and 
'pplicabiJ ity  I o elastic  shells  and  fluid  sheets. 


'•'«  Sec  nn 
® fi  Section 


Tirr  - 


out.  by  F.  and  F.  I'oscernl  I 41.  In  their  work,  the  rosceruts  represented  the 
orientation  of  each  point  of  their  continuum  by  a set  of  mutually  perpendicular 
rigid  vectors.  The  purely  kinematical  aspects  of  oriented  bodies  characterized 
by  ordinary  displacement  and  the  independent  deformation  of  N deformable  vectors 
in  N-dimensional  space  has  been  discussed  by  Kricksen  and  Truesdell  t V>  1 » who  also 
introduced  the  terminology  of  directors . A complete  general  theory  of  a Cosserat 
surface  with  a single  deformable  director  given  by  Green,  Naghdi  and  Wainwright 
[b]  was  developed  within  the  framework  of  thennomochanics ; and  their  derivation 
in  to)  is  carried  out  mainly  from  an  appropriate  energy  equation,  together  with 
invariance  requirements  under  superposed  rigid  body  motions.  A related  develop- 
ment utilizing  three  directors  at  each  point  of  the  surface,  in  the  context  of 
a purely  mechanical  theory  and  with  the  use  of  a virtual  work  principle,  is  given 
by  Cohen  and  DoSilva  [7].  A further  development  of  the  basic  theory  of  a Cosserat 
surface  is  given  by  Naghdi  f 1 1 , which  also  contains  additional  historical  remarks 
relevant  to  oriented  continua  and  th<  development  of  the  theory  of  elastic  shells. 

Before  describing  the  contents  of  this  lecture,  it  is  helpful  to  provide  a 
definition  for  a shell-like  body.  To  this  end,  consider  a three-dimensional  body 
embedded  in  a Euclidean  "2-space  and  let  its  particles  be  identified  by  a con- 


vected  coordinate  system^  01,  (i = 1,2,3).  Let 

p = p(ex,e2,e’,t) 


(1.1) 


denote  the  position  vector,  relative  to  a fixed  origin,  of  a typical  particle  of 
the  three-dimensional  body  in  the  present  configuration  at  time  t.  In  a 
reference  configuration  of  the  body,  not  necessarily  the  initial  configuration, 
we  denote  the  position  vector  by 


P 


Ue1  ,e2,e3) 


(1.2) 


If  the  reference  configuration  of  the  three-dimensional  body  is  specified  by  the 
initial  configuration,  ry  at  time  t-0,  then  the  right-hand  side  of  (1.2)  can 


be  identified  with  p(9l,0). 
notation 


Now,  for  convenience,  set  Q-*’  = § and  adopt  the 


e1  = (eV,) 


(1.3) 


K*  t-ping  the  above  in  mind,  we  begin  by  defining  first  in  descriptive  terms  what 
is  meant  by  1 thin  shell  or  a shell-like  body  in  its  reference  configuration. 
Consider  a two-dimensional  surface,  called  a reference  surface,  which  may  be 
defined  by  the  parametric  equation  § = 0;  and  let  R and  A^  denote,  respectively, 
the  position  vector  and  the  unit  normal  to  the  reference  surface.  Imagine  now 
material  filaments  from  above  and  below  surrounding  the  surface  along  the  normal 
at  each  point  of  the  reference  surface.  Suppose  further  that  the  bounding 
surfaces  formed  by  the  end  points  of  the  material  filaments  are  equidistant  from 
the  reference  surface.  Such  a three-dimensional  body  is  called  a shell  if  the 
dimension  of  the  body  along  the  normals,  called  the  thickness  and  denoted  by  h, 
is  cmal 1 . A shell  is  said  to  be  thin  if  its  thickness  is  much  smaller  than  a 
certain  characteristic  length  of  the  reference  surface  such  as  the  minimum 
radius  of  curvnt'ire  of  the  reference  surface.*  If  h is  constant,  the  shell  is 
said  to  be  of  uniform  thickness,  otherwise  of  variable  thickness.  A three- 
dimensional  shell-like  body  in  its  reference  configuration  may  be  depicted  as  in 
Fig.  1,  where  % * 0 is  shown  as  the  middle  surface  and  the  bounding  surfaces  are 
specified  by  5 - - h/2  with  h a constant.  Let  the  region  of  space  occupied  by 
the  shell  in  tin  reference  configuration  bo  covered  by  n norma 1 coordinate 
r.yri  cm  (9°\£)  with  Q along  the  normal  to  the  middle  surface  and  0a  on  the 
r*dor>  nce  surface  (£-£-0).  Then,  the  position  £ of  any  point  of  the  shell  in 


Recall  that  when  the  particles  of  a continuum  are  referred  to  a convected  co- 
ordinate system,  the  numerical  values  of  the  coordinates  associated  with  each 
particle  remain  the  same  for  all  time. 


^Tn  the  case  of  a plate,  since  the  reference  surface  is  a plane,  the  character- 
istic length  is  taken  to  bo  the  smallest  dimension  of  the  reference  plane. 


2. 


its  reference  configuration  may  be  expressed  as  a function  of*  (9a'»C): 

£ = = K(e“)  + £A3(e“)  . 

With  the  above  description  in  mind  we  now  proceed  to  provide  a more  genera i 
definition  for  shell-like  bodies  in  fairly  precise  terms.  Since  a material 
surface  in  the  three-dimensional  body  can  be  defined  by  the  equation  £ = £(0a)> 
it  follows  that  the  equations  resulting  from  (1.1)  and  (1.2)  with  5=5(0*) 
represent  the  parametric  forms  of  the  material  surface  in  the  present  and  the 
reference  configurations,  respectively.  In  particular,  the  equation  §-0  define: 
a surface  in  space  at  time  t,  which  we  assume  to  be  smooth  and  nonintersecting. 
Any  point  of  this  surface  is  specified  by  the  position  vector  r , relative  to  the 
same  fixed  origin  to  which  p is  referred,  where 

r = r(e‘  ,9'')  = p(eJ  .e^.o.t)  . (1.5) 

Let  the  boundary  of  the  three-dimensional  continuum  be  specified  by  the  material 
sur  faces 

? = ?2(eV)  , i = 51(e1.e2)  > §2  < g1  , U-6) 

with  the  surface  * = 0 lying  entirely  between  them,  and  a material  surface 

ffeLe2)  = o , (1.7) 

which  is  chosen  such  that  5 = const,  form  closed  smooth  curves  on  the  surface 

(1.7) . 

by  way  of  additional  background  information,  suppose  now  that  £ in  (1.1  ) is 
a continuous  function  of  Q^t,  and  has  continuous  space  derivatives  of  order  1 
and  continuous  time  derivatives  of  order  2 in  the  bounded  region  ^ • 

Hence,  to  any  required  degree  of  approximation,  £ bo  represented  as  a 
polynomial  in  - with  coefficients  which  are  continuous  functions  of  0^,t. 

Hcwvvr , instead  of  considering  a general  representation  of  this  kind,  we 
restrict  attention  here  to  the  approximation 

p = r*§d,  (1-' 

wh- r-  £ is  defined  by  (1.‘  and  d d(0*,t).  The  expression  corresponding  to 

(1.8)  in  the  reference  configuration  is 

P = F + §£  , (1.9) 

where  £ £(  0°  j is  the  reference  value  of  (1.5)  and  D = D(flP)  is  the  reference 
value  of  i.  fince  a representation  of  the  form  (l .4)  can  always  b*  chosen  in 
any  one  configuration  without  loss  in  generality,  it  is  of  interest  to  indicate 
the  relationship  between  the  right-hand  sides  of  (1.9)  and  (l.U),  let  D be 
specified  long  the  normal  to  the  surface  5 = 0,  i.e.,  let 

D = DA^  , (1.10) 

where  l1  is  the  magnitude  of  £.  Then,  after  equating  (l.M  and  (1*9)>  we  obtain 

C - (1.11) 

tir  t r a ns.  format  ion  r«*iatii-n  b-twoen  £ arid  5.  Hence  for  a she!  1-like  body, 
w i . I ■ ‘ ( ! . 8 ) represents  an  approximation  for  the  position  vetor  in  the  current 
configuration  at  time  t,  the  simple  representation  (1.9)  in  the  reference  con- 
figuration with  I)  specified  in  the  form  (1.10)  involves  no  loss  in  generality. 

In  the  rest  of  the  paper,  we  consider  a fairly  general  development  of  tin* 
theory  oi  shells  by  direct  approach  suitable  for  application  to  shell-like 
bodies  as  defined  in  the  preceding  two  paragraphs.  Specifically,  the  next  two 


For  details  see  [1,  pages  kU2  and  1*71]. 


sections  deal,  respectively,  with  kinematics  and  the  baric  principles  for  shells 
based  on  a Cosserat  surface  and  section  4 is  concerned  with  a discussion  of 
constitutive  equations  for  elastic  shells.  The  final  section  of  Ui»*  paper  con- 
tains a brief  siamnnry  of  some  relevant  researches  pertaining  to  elastic  shells. 
An  effort  is  made  to  empioy  a direct  (coordinate  free)  notation  in  keeping  with 
recent  trends  in  finite  elasticity  and  the  subject  of  this  symposium.  r Otherwise 
our  notation,  with  some  minor  differences,  is  essentially  that  in  fl).^  Through 
out  the  paper,  Latin  indices  (subscripts  or  superscripts)  have  the  range  1,2,3, 
Greek  indices  have  the  range  1,2  and  the  usual  summation  convention  is  employed. 

2.  KINEMATICS 


Deformable  media  which  are  modelled  by  a material  surface,  embedded  in  a 
Euclidean  3-space,  together  with  K (K = 1,2 , . . . ,N)  deformable  vector  fields  -- 
called  directors  --  attached  to  every  point  of  the  material  surface  are  called 
Coes or at  surfaces  or  directed  surfaces  and  may  be  conveniently  referred  to  as  cK 
In  the  absence  of  the  directors,  we  merely  have  a two-dimensional  material 
surface  which  can  serve  as  a model  for  the  construction  by  direct  approach  of 
the  membrane  theory  of  shells.  With  K = l,  the  directed  medium  is  a body  Cj_  C 
consisting  of  a material  surface  and  a single  deformable  director  attached  to 
every  point  of  the  material  surface  of  C*  The  latter  is  the  simplest  model  for 
the  construction  of  a general  bending  theory  of  thin  shells;  and,  for  simplicity 
we  restrict  attention  to  this  particular  model  in  the  discussion  that  follows. 6 
L-t  the  particles  of  the  material  surface  of  C be  identified  by  a system  of 
converted  coordinator,  Q01  (a-  1,2);  let  the  material  surface  of  C in  the 

present  configuration  at  time  t,  be  described  by  its  position  vector  r relative 
to  a fixed  origin;  let  and  denote,  respectively,  the  base  vectors  along 
the  B'^-curves  on  J and  the  unit  normal  to  and  let  d stand  for  the  director  at 
£.  Then,  a motion  of  the  Cosserat  surface  is  defined  by  vector-valued  functions 
which  assign  position  r and  director  d to  each  particle  of  C at  each  instant  of 
time,  i.e.,' 


r=r(ea,t)  , d=d(9a,t) 


[a^a^d]  / 0 


(2.1) 


whore 


a (0\t) 
~or 


dr 


d9 


(2.2) 


and  the  condition  12-1 ) ^ ensures  that  the  director  d is  nowhere  tangent  tc  J. 

It  is  convenient  to  introduce  here  a slightly  different  notation  than  that 
adopted  in  Naghdi  [1]  and  a number  of  previous  papers.  Thus , we  put 


d = a 

—a  ~a 


i3  = £ 


(2. 


and  observe  that,  in  view  of  (2.l)j  and  (2. 3),  d^,dg,dg  ore  linearly  independent 
vectors.  Hence,  we  may  introduce  a set  of  reciprocal  vectors  d1  such  that 

d.  • dJ  = 6J’  (2.4  ' 

where  b'l  is  the  Kronecker  symbol  in  three-space.  Whenever  desirable,  the  nota- 
tions L - (d ^ ,<L> ,d^ ) and  a^,d)  will  be  used  interchangeably  throughout  the 
paper  depending  on J the  particular  context. 

Consider  now  a reference  configuration,  not  necessarily  the  initial 


; I’ll* • coordinate  free  notation  employed  here  is  similar  to  that  used  by  Carroll 
and  Naghdi  [8]  but  some  of  the  definitions  are  different. 

For  a more  general  development  of  the  theory  of  Cosserat  surfaces  see  for 
example  section  2 of  a recent  paper  by  Green  and  Naghdi  [9]  which,  however, 

!•  als  with  fluid  sheets  and  water  waves. 

7 

For  convenience,  we  adopt  the  notation  for  r in  (1.5)  and  (1.8)  also  for  the 
surfac*  '2.l)j.  This  permits  an  easy  identification  of  the  two  surfaces,  if 
desired. 


configuration,  of  the  Oosserut  surface  C»  In  the  reference  configuration,  let 
the  material  surface  of  C he  referred  to  by  $ with  K as  its  position  vector;  let 
A^Aj  denote,  respectively,  the  base  vectors  along  the  ©^-curves  on  '»  and  the 
unit  normal  to  *V,  and  let  D be  the  director  at  K.  Then  corresponding  to  (2.J), 
we  have 

R = R(8a)  , £=£(ea)  , ^ o , (?.',) 


A - A (9V)  = (2.6) 

^ df 

and  (2.^)^  ensures  that  P is  nowhere  tangent  to  the  rurface  If  the  reference 
configuration  of  C is  specified  to  be  the  initial  configure  * ion,  ray  at  time 
t o,  then  the  vector-valued  functions  on  the  right-hand  sides  ol  (2.^)i  . can  b* 

identified  with  r(0a,>)  and  d(0o,O),  respectively.  Analogously  to  (2.3;  J*' wo  set 

1>  = A , P - D (2.r) 

~q  ~a  ~3  ~ 


and  note  that  the  dual  of  (2.4)  is  given  by 


As  in  [8],  we  introduce  the  notations  grad  and  Grad  to  denote  the  right 
spatial  and  material  gradient  operators,  respectively,  with  respect  to  the 
position  on  the  surface  J in  the  current  configuration  and  on  the  surface  **  in 
the  reference  configuration.  The  corresponding  divergence  operators  will  be 
denoted  by  liv  • nd  Div,  respectively.  In  particular,  for  a vector-valued  func- 
tion V^fjPjt),  we  have 


grad  V V &a 

- ~,a  ~ 

Grad  V V 8 A0 

~ O'  ~ 


div  V = V 


Div  V - V • A 

~ ~,a  ~ 


>«e  d fill'*  measure  of  deformation  by  the  tensor  F,  n:uneiy 

F .1.8 1)1  - Grad  r + d.gD1  , (2.10) 

~ ~i  ~ — ■ ~3  ~ 

where  the  symbol  S denotes  the  tensor  product.  Keeping  the  notations  (2-3)  and 
(2.7 ) in  mind,  we  observe  that 

F D = F A = a = d , 

~~ a ~Q  -or  (2.11) 

F£^  = FD==d  = d^  . 

i.v  tne  definition  of  the  determinant  of  a second  order  tensor  and  the  conditions 
(2.1  i 4 < nd  (2.5) y as  well  as  the  fact  that  for  continuous  motions  the  scalar 
triple  products  in  (2.1)  and  (2. 5)  must  have  the  same  sign,  we  obtain 

det  ¥ = [d^dJ/tDJ^Dg]  > 0 . (2.12) 

Ti  e tensor  F,  a linear  operator  on  vectors  in  3-space,  is  nonsingular ; and  tliere 
exists,  therefore,  the  inverse  deformation  gradient  tensor  F“^  defined  by 


Tii*'  inverse  operator  £“J  transforms  vectors  in  the  present  configuration  into 
vectors  in  the  reference  configuration,  i.e., 

The  symbol  F in  the  paper  of  Carroll  anil  Nnghdi  f8]  stands  for  a different, 
quantity.  The  term  Grad  r in  (2.10)2  corresponds  to  the  deformation  gradient 
tensor  F in  [8], 


(2-Ut  i 


f‘1d.  - I). 

~ ~»1  ~1 

and  it  foil owe  that 

f*1K  = F F"1  --  I = <1.  Sd1  = l>.  SiD*  , (2.1'.) 

~ ~ ~ ~ ~ ~1  ~ ->-l  ~ 

where  £ is  the  unit  tensor  in  3-space.  We  also  introduce  here  the  director 
gradient  tensor  by 

G = Grad  d = d„  8D°  = d »Da  . (2.16) 

~ ~ ~-3 , ex  ~ a ~ 

The  velocity  and  the  director  velocity  vectors  at  a point  of  J and  at  time  t 
are  defined  by 

v_=r(o“,t)  , w=d(9“,t)  , (2.17) 


where  a superposed  dot  denotes  tin*  material  time  derivatives  with  respect  to  t 

holding  0**  fixed,  fince  a - v , we  have 

~,0' 


F = d.  8D1  = d ®Da  + d *D3  = v 8Da  + w®D3 
~ ~ ~o  ~ ~3  ~ ~,Q  ~ ~ ~ 

G = d,  »1)“  = w *1)®  . 

~ ~3,  a ~ ~.a  ~ 


(2.18) 


Also, 


F V ^ = d^d1  = grad  v 

* -l  & 

I F = w $$d  grad  v 
~ ~ ~»a  ~ ~ 


(2.19) 


3.  BASH’  HilNt’IPLfiS 

In  the  development  of  this  section,  we  follow  the  mode  of  derivation  of  the 
basic  theory  employed  by  Naghdi  1,  Section  ft].  let  P,  bounded  by  a closed 
curve  dP,  be  -my  part  of  J occupied  by  an  arbitrary  material  region  of  the 
surface  of  r in  the  present  configuration  at  time  t and  let 

v = v“d  S V d“  (3.1) 

~ ~a  a~ 

bo  the  outward  unit  normal  to  9p.  Wo  now  introduce  various  field  quantities 
issocia' : with  tie  Coss.-r-tt  surface  as  follows:  The  mass  density  p-  p(0^,t)  of 
t.li"  surfae.  J in  tii'.'  present  configuration;  tlr  internal  mechanical  action  across 
any  curve  dp  characterised  by  the  contact  fore" * n~n(9a,t:v)  and  the  contact 
director  force  m mto^jt  ; v , each  measured  per  ’ mif~l earth  of  a curve  in  the 
present,  configuration;  the~intrinsic  (surface)  director  force  k = k(9a,t)  per  unit 
art  a of  J.  which  makes  no  contribution  to  the  supply  of  moira  nt~o f "momentum ; the 
issigned  force  £ f(9'),t  ' and  the  assigned  director  force  l £(9a,t),  each  jar 

unit  mass  of  th"  surface  J\  and  th‘  inertia  coefficients  ajq  <»|ij(9')>  (M-1,2), 
are  assumed  to  be  independent  of  time.  We  assume  that  the  kinetic  energy  or  the 
'oss- rat  surface  C per  unit  area  of  J in  the  present  configuration  is  given  by 

T - 1 p(v  . v + 2ttiV  • w + a2w  . w)  . (3.2) 

'Wv  further  define  the  momentum  corresponding  to  the  velocity  v and  the  director 
momentum  corresponding  to  the  director  velocity  w by 

p(X.tQ,j”)  > § = P(“g+0^iL)  • (3-3) 


Hi.  notations  for  the  contact  force  n,  the  contact  director  force  m,  and  the 
r'irface  director  force  k differ  fron~Naghdi  [lj  and  the  previous  papers  on 
tf.*  subject.  In  fact,  the  vector  fields  n,.m,k  of  the  present  paper  correspond, 
respectively  to  N_,M,m  in  Naghdi  (!]. 


•>. 


L 


= r p(r  x f + d x £)da  + [ (r  x n + d x m)ds  , 
JP  ~ ~ Jdf>~  ~ ~ ~ 


wh* re  da  is  the  element  of  area  and  ds  is  the  line  element  in  the  surface  J. 

Th*  first  of  (3.5)  is  a mathematical  statement  of  the  conservation  of  mass, 
th*  second  that  of  the  linear  momentum,  the  third  that  of  the  director  momentum 
and  the  fourth  is  the  conservation  of  moment  of  momentum.  The  basic  structure 
of  (3.5)2  **nd  (3.5)o  and  their  forms  are  analogous  to  the  corresponding  conserva- 
tion laws  of  the  classical  three-dimensional  continuum  field  theory.  Th** 
structures  of  (3.5)0  and  arc  less  obvious,  but  a motivation  for  their 

forms  is  provided  by  a derivation  of  the  basic  field  equations  for  shell-like 
bodies  obtained  from  the  three-dimensional  equations  of  continuum  mechanics  in 
which  the  position  vector  £ in  3-space  is  approximated  by  an  expression  of  the 
form  ( 1 . 8 ) . 

It  should  be  noted  her*  that  the  conservation  laws  (3.5)  ar**  consistent 
with  the  invariance  conditions,  under  superposed  rigid  body  motions,  which 
ordinarily  have  wide  acceptance  in  continuum  mechanics.  Moreover,  as  shown  in 
. S :tion  J , th  eon.  rvat  ion  lawj  • . • , • . an  1 . , equi  Lent 

to.  and  c**n  be  derived  from,  an  appropriate  conservation  of  * n*  r *y  and  t in- 
variance conditions  under  superposed  rigid  body  motions.  Hie  conservation  law 
(3.5)'  for  the  director  momentum  must  be  postulated  separately. 

Returning  to  the  conservation  laws  (3*5)  > we  note  that  under  suit  a 1 • 
continuity  assumptions  the  contact  force  n and  the  contact  director  force  m can 
be  expressed  in  the  forms  (for  details  see  i 1 , Section  8jv): 


Depending  on  the  choice  of  the  physical  dimension  of  d and  with  reference  to 
m ,1  and  k,  the  terminologies  of  the  contact  director  couple,  the  assigned 
iirector  couple  and  the  intrinsic  director  couple,  respectively,  are  also  used 
in  the  literatur*  . In  particular,  the  latter  terminologies  are  employed  in 
ill,  where  d is  taken  to  be  dimensionless. 

1 ‘As  the  integrals  on  the  left-hand  sides  of  (3.5)^  3 ^ allow  for  coupling  in 
inertia  terms,  they  are  slightly  more  general  than  £hc  corresponding  expres- 
sions in  1 1 1 . The  conservation  laws  (3»5)  with  the  coefficients  - 0, 

02s  reduce  to  those  given  by  equations  (8.17)  in  [1]. 


■r 

V ft  V 
a 


n N^v  N*v  y in  * M°v  = M v > (3.6) 

where  the  second  order  tensors  N,M  and  the  vectors  are  related  through 


N - d »N 
~ ~a  ~ 


N“  = HTdB 


M = d ®M“  , = MTd“ 


(3-7) 


and  the  superscript  T denotes  transpose.  Also,  for  convenience,  we  introduce  a 
tensor  v tiirough 


j i = 


K = KTd3 


(3-8) 


V.'Lth  th**  use  of  (3*6)  and  by  usual  procedures,  from  the  conservation  laws  (3.5) 
follow  the  local  equations  in  either  of  two  equivalent  forms: 


pa*  = X or  p + pa  • v =0 
~ ~>C* 

+ p£  = (3  ( + “l^.)  > 


pi-k  - p(a  v + or0w ) , 


a x Na  + d x k + d x M01  = 0 
~o  ~ ~ ~ ~,ar  ~ ~ 


(3-9) 


P * p 


div  v = 0 


div  H*  t pi'  = Ptv^+a^w)  , 
div  K1  i o!-k  --  p(a,  v e ct,w) 

(N  + K + G r'1  M)  = (N+  K + G F'1  M)1 


(3-10) 


win  re  \ is  n function  of  0 only,  a conroa  denotes  partial  differentiation  with 
respect  to  0"*,  a vertical  line  stands  for  covariant  differentiation  with  respect 
to  the  metric  tensor  of  the  surface  J and 

a?  = la  a a J . ( 1.11 ) 

A . u,  i'..  tile  definition  of  tut*  rirnt  di . --ivnc*'  ui'  *t  tensor  field,  wv  have 

iiv  \ • Uv  if  if |ff  • 

II  ir  interesting  that  the  last  statement  in  (3.1  1 ' is  similar  to  the  symmetry  of 
'li  - stress  in  the  Hire,  -dimensional  theory.  Furthermore,  it  may  be  observed  that 
i*'Z.  ani  d a x If*  ar-,  respectively,  the  axial  vectors  of  l K1  - j 

-ni  i/i  f’/i.  r:,;  hh:\:  ' ~ 

At  £nis  point,  ir  is  desirable  to  make  some  observations  concerning  the  two 
s-  t"  of  field  equations  (3.  “)  and  (3.10).  If  one  assumes  an  expression  of  the 
form  1.  for  the  position  vector  of  the  sh-’ll-like  body,  then  field  equations 
of  the  same  form  as  (3.9'  can  also  be  derived  from  the  ttiree-dimensional  field 
equations  of  classical  continuum  mechanics  by  suitable  integration  with  respect 
to  x between  the  limits  and  ?2  defined  in  (1.6).  In  a derivation  of  this 
hind,  one  also  needs  to  define  appropriate  resultants  corresponding  to 
as  well  as  certain  integrals  corresponding  to  the  Inertia  coefficients .~ 

I he  assigned  field  £,  which  occurs  in  the  equations  of  motion  (3.0)^  and 
> 3.1  • • represents  the  combined  effect  of  (i)  the  stress  vector  on  the  major 
surfaces  of  the  shell-like  body  denoted  by  f and  (ii)  a contribution  arising 
from  the  three-dimensional  body  force  denoted  by  f.  . A parallel  statement  holds 
for  the  assigned  field  £ in  (3-9)1*  and  (3.10)j.  Thus,  we  may  write 


f = f,  + f 


i ~ L + l 


(3-13) 


llu*  various  quant, itios  in  (3.1:0  arc  lVeo  to  he*  specified  in  a niMirur  which  de- 
pends on  the  particular  application  in  mind  and,  in  the  context,  ol'  tin*  theory  of 
a i‘osserat,  surface,  the  inertia  coefficients  require  constitutive  equations. 
Indeed,  both  f^,  and  , as  well  as  f^  and  can  be  identified  with  correspond- 
ing expressions  in  a derivation  from  the  three-dimensional  equations  (for  details 
see  [1,2]).  Likewise,  the  inertia  coefficients  may  be  identified  with  easily 
accessible  results  from  the  three-dimensional  theory. 

Both  sets  of  equations  (3-9)  and  (3.10)  are  simple  in  appearance  but  they 
conceal  the  relative  complexity  of  the  results.  Although  it  is  the  component 
forms  of  the  equations  of  motion  that  are  useful  for  application  to  specific 
problems,  on  occasions  the  forms  (3.9)  and  (3.U>)  offer  son-  advantages, 
especially  in  discussions  pertaining  to  aspects  of  the  general,  theory,  kither 
of  the  two  sets  can  be  expressed  in  terms  of  tensor  components  referred  to  either 
the  basis  a.  - (a_,a^ ) or  the  basis  d-  (<k  ,d^)  defined  by  (2.3).  We  do  not  dis- 
play result^  of  this  kind  her**  but  note  that  th**  component  forms  of  (3.*))  in 
terns  of  tensor  components  referred  to  a^  are  presented  and  utilized  in  [lj. 

Before  closing  this  section,  w«>  recall  that  the  rate  of  work  by  ail  contact 
and  assigned  forces  acting  on  the  part  p and  on  its  boundary  dP  minus  the  rate 
of  increase  of  the  kinetic  energy  in  p can  be  reduced  to 


where 


J p(f  • v + 1 • w )da  + J (n  • v + m . w 


)ds 


dP 


d 1 r 

” dt  J P ( v ‘ v + 2ar^v  • w + OgW  • w)d.a  = J 


P do 


(3-11*) 


V - • v + k • w + Ma  • v 

~ ~ ~ ~ ~,a 

tr{:;  ('.rad  v t k(w  ®it’)  t M (.'rad  w] 

- tr{[  (k+k)k  + m ajK-1]  . 


(3.15) 


I*.  KI.AP.TIC  SHELLS 


Within  the  scope  of  the  general  theory  of  a Cosserat  surface,  we  discuss 
briefly  the  constitutive  equations  for  elastic  shells  in  the  presence  of  finite 
lei  rmoti  n.  reliminary  to  the  discussion  that  follows,  we  assume  the  existence 
of  a strain  energy  or  stored  energy  per  unit  mass  f=*(0  ,t)  such  that  pi  is 
equal  to  the  mechanical  power  defined  by  (3*15),  i.e., 

P = p^f  (,  U . l 


In  passing,  we  note  that  by  us*-  of  (U.l)  and  the  definition  of  the  strain  energy 
fur  a foss-rat  surface  C,  namely 


U = f p ^ da  , (*1.2) 

JP 

f ra  m (3*1 0 follows  the  analogue  of  a well-known  result  in  the  three-dimensional 
theory:  th*-*  rate  of  work  by  the  contact  and  the  assigned  forces  and  director 
forcer  acting  on  P and  on  its  boundary  dP  is  equal  to  the  sum  of  the  rate  of 
kinetic  energy  in  P and  the  rate  of  the  strain  energy  in  p. 

Be turning  to  our  main  objective  in  this  section  concerning  the  development 
ot  nonlinear  constitutive  equations  for  elastic  shells,  we  assume  that  the 
strain  energy  density  ifc  at  each  material  point  of  C and  for  all  t is  specified 
by  a response  function  which  depends  on  £ ,d  and  their  partial  derivatives  with 
r-  ; n*  ct  to  Qa . But  since  the  response  function  must  remain  unaltered  under 
up*. -i*p< • I ri  id  body  translational  displacement,  the  dependence  on  r must,  b*- 
* x(*Ji..iei.  Ihuc,  the  conr titutive  assumption  for  the  strain  onerg y density  can 


be  written  as 


* 


* '(r  ,d,i!  ; X) 

~>a >£» 


•:v.i  w*'  also  make  similar  constitutive  assumptions  l*or  in  these  con- 

stitutive equations,  which  represent  the  mechanical  response  of  the  medium,  the 
dependence  of  the  response  functions  on  the  local  geometrical  properties  of  a 
reference  state  and  material  inhomogeneity  is  indicated  through  the  argument  X. 

A general  development  of  various  aspects  of  constitutive  theory  of  elastic 
shells  based  on  assumptions  of  the  type  (4.3)  °r  variants  thereof  is  given  in 
[l,  Section  13).  In  the  rest  of  this  section,  we  limit  the  discussion  to  an 
elastic  shell  which  is  homogeneous  in  its  reference  configuration  and  suppose 
also  that  the  dependence  of  the  response  functions  on  the  properties  of  the 
reference  state  occurs  tlirough  the  values  of  the  kinematical  variables  in  the 
reference  state.  Then,  in  place  of  (4.3),  we  have 


V = 


r »d,d  „ 
a , at 


A , 

~o 


D,D  ) 

, a 


(4.4) 


with  similar  assumptions  for  N^k^M^.  If  instead  of  the  kinematic  variables  n 
(4.4)  the  response  functions  are  assumed  to  depend  on  Grad  r^,d,Grad  d and  their 
reference  values,  then  we  may  write 


■]i  = ^(Gr;i  r,d,Grad  D ; Grad  h,D,Grad  D) 


•• 


with  similar  assumptions  for  N ,k ,M.  Still  another  form  for  the  response  func- 
tions could  involve  the  tensors  V,G.  In  this  connection,  we  observe  that  the 
kinematic  measure  F in  (2.10)  alreadyjinvol ves  the  reference  values  and  h :o 
that  a response  function  of  the  form  ^ may  be  expressed  as  a different  function 
of  F,  I and  the  reference  value  of  Thus,  corresponding  to  the  assumption  (4.4  . 
we  may  write 


* = ?(£>£ : 

Wh'-Te 


(4.6) 


„G  = Grad  D = D 8D®  , (4.7) 
along  with  similar  assumptions  for 

Keeping  the  above  constitutive  assumptions  in  mind,  with  the  use  of  (4.1) 
nd  (3.15)i  p 3,  by  usual  techniques  we  obtain  the  following  alternative  forms 
for  the  constitutive  equations: 


or 


or 


£ ' p la  ’ 

m“  = P -St- 

~ M dd 

(4.8) 

M )T 

d Grad  r 

> £ p |j 

, M - - p Grad 

(4.9) 

P(  i X f*  ; 

-Of  ~ o. 


K p(d  8 


D3)(§)T  , M = oF(  — 


^ 0 Grad  r (“4  1 . (t . . 

dG  ~ o'.  1 


In  addition,  the  strain  energy  response  function  in  (4.8)  is  restricted  by  the 
invariance  condition1-^ 


x + d x + d x ■ = 0 

,or  dr  ~ dd  ~,  cr  dd  ~ 

~yOt  ~ 


while  the  corresponding  restriction  for  ^ is 

£ - £'  > £ - ^ (f  )'r  ♦ (dg  * £3  ) (f  >T  ♦ G(f  : 


(4.11). 


(4.12) 


12 

‘This  form  was  employed  in  the  paper  of  Carroll  and  Nnghdi  [8]. 
J Compare  the  combination  of  (4.8)  and  (4.11)  with  ( 3 • 9)|^ - 


1 . 


A condition  similar  to  (4.  I,’ )j  holds  for  t.li-  function  ^ in  ,4.  **  . 

V.e  do  not  discuss  here  the  reduced  forms  of  the  above  eonstituli v-*  equations 
resulting  from  invariance  requirements  under  superposed  rigid  body  motions,  but 
for  such  reductions  refer  the  reader  to  [1,  Section  13].  Just  as  with  th« 
equations  of  motion,  it  is  necessary  in  applications  to  specific  problems  to 
obtain  alternative  forms  of  the  above  constitutive  equations  or  their  reduced 
forms  in  terms  of  tensor  components.  Such  component  forms  may  b<-  express*  .1  with 
respect  to  bases  n,j  or  d • , or  corresponding  bases  in  the  reiVren-'e  confi»?ira- 
tic-n.  Reduced  for  u f (4.d)  have  been  utilized  extensively  in  I 1 , Chapters  D 
and  Kj . 

. ADDITIONAL  REMARKS 

In  this  section  we  briefly  curanent  on  some  : pec  ini  ct..*s  * f the  - n.  ra  L 
theory  and  also  mention  some  recent  research'  s which  Lear  on  th--  .arious  aspeetr 
' istic  shells . Although  these  ievelopmenti  will  • *ri  t ntext 

ot  a mechanical  theory,  some  of  the  reference's  cited  are  : -t:.*  r * mo r--  g-neral 
and  include  thermal  effects. 

As  noted  in  Section  i,  the  well-known  ra  s.i  t r . • > i p,  : . •.  ■ : • 

; ■ ined  is  a s\  ec i a 1 case  >f  t he  ner  tl  rj  * 

suppressing  the  effect  of  the  director  and  curr-spend ing  kin--?  • variabl- -s  and 

this  is  discussed  briefly  in  ; 1 , Section  i4] . A levelopment  1 nothei  special 
theory,  known  as  the  i nex t < ns  i < >na  1 thet  ry , w erein  t Length 
th  surface  of  C is  assume  1 t o rema in  constant  thi  igh  it  i • . • . . i 

^ ■ i in  [l.  Section  L4] . Sii  Llarly,  a non  ar  i str t 1 tl  * 

by  direct  approach,  motivated  mainly  by  the  classical  tin' cry  corr  sponding  to 

Kirchhoff-Love  theory  of  shells  and  plates,  is  r iven  by  Naga  li  f'-ctions  1. 
s i 15].  & constrained  th  ry  of  an  elastic  rat  surl  :ussed  b> 

ire- n and  L’aghdi  [10]  and  includes  as  a special  cas-*  a theory  which  is  in  1-1 
correspondence  with  the  restricted  theory  mentioned  above. 

The  nonlinear  constitutive  equations  in  Lection  4 are  v.ali  i for  an  elastic 
Cos rerat  surface  which  may  L-  anisotropic  with  refer- nee  to  preferred  directions 
assoc: aled  with  material  points  of  C-  A general  discussion  of  mat-rial  symme- 
tries for  ch<  Lis  is  . , Seel  n 13].  ' rr<  . id  N hdi  ha « 

subsequently  ex'imined  tii--  influ*'nee  of  the  r-'fcr-nce  g.eometry  on  the  response  of 
•‘.•■stic  shells  by  assuming  the  ••xist*,nce  of  a local  preferred  state  of  the  body 
' : ; ■ ting  t ha f t . ! i 1 the  reference  geometry,  as  in  (4.4), 

■ *ears  through  the  va  I .•  ; : • as.s  t i tut  iv*  variables  in  the  pr- for  red  state. 

Matei  resti  t • stic  shells  ha\  1 n discussed  Lso,  fron 

• different  point  of  view,  : '-r:  n iLlj  who  has  indicated  in  [.12]  a compari- 

son with  the  results  contain-  i :r. 

oner  f w • I i t hells,  based  on  th< 

' ■ ' ' • ■ n discussed  1 iricksen  . A related 

study  on  ti  • subject,  limited  ...  * ■ w propagation  in  a surface  not  endowed 
with  ■ dir-ctor,  was  giv-  r:  • arl’-r  : * ?.•  s and  Suh  [l4]. 

• ry  of  su  • ■ r posed  Larg  formation  of  at 

elastic  * • ■ rat  surface , -iu*  - r-  at-  i problems,  of  stability  and  vibrations 
1 ' rn  l plates,  1 nd  Naghdl  . : • 

system  of  linear  equations  *h  s a •*  • ris  i ng  th*-  initial  mixed  boundary-value 
prct  .em  of  elastic  shells,  haghdi  i I r * : ; -•!  have  obtained  a unique n-u  s 

without  th  . f definitenei  for  the  strain  en<  rgy  :•  ns i ty . 

Th'  theorem  in  loj  holds  for  nonhomcg-*neous  and  anisotropic  shells  undergoing, 
small  motions  superposed  on  a large  deformation. 

In  still  another  study,  th-  theory  of  a 'ossernt  surface  has  been  employed 
»;>  ffhdi  [17]  t.o  formulate  contact  problenu  : nd  tes.  In  the 

derivfition  of  sh-'ll  theory  from  the  three-dimensional  equations,  equations  of 
motion  in  terms  of  resultants  and  detailed  consideration  of  constitutive  equa- 
tions for  shells  are  usually  obtained  relative  to  an  interior  surface,  rather 
than  one  of  the  major  surfaces  of  the  shell -like  body  which  may  be  the  contacting 
surface;  the  interior  surface  ordinarily  is  identified  with  the  middle  surface  of 
th-  rh.  il  or  plate  in  the  reference  configuration.  In  the  development,  of  shell 
theory  by  direct  approach,  although  the  material  surface  of  C may  be  id.-ntifiod 
with  any  surface  of  the  (three-dimensional)  rhell-like  body,  nevertheless  the 


1 1. 


complete  discussion  of  constitutive  equations  and  the  identi fication  of  the 
inertia  coefficients  and  the  assigned  fields  again  require  explicit  use  of  a 
reference  surface  in  the  shell-like  body.  For  certain  problems  it  is  more  natural 
and  conceptually  more  appealing  to  select  one  of  the  two  major  surfaces  as  the 
reference  surface  but  then  the  detailed  available  development  of  the  constitutive 
equations,  as  well  as  identification  of  such  quantities  as  the  inertia  coeffi- 
cients, have  to  be  reconsidered  relative  to  the  new  surface.  This  problcsn  is 
resolved  in  [17]  by  deriving  appropriate  transformation  relations  which  relate 
th'-  kinetic  variables  N^jk,^  (and  hence  the  response  functions  ) in  the  two 
formulations.  The  results  in  [17]  are  applicable  to  any  shell -like  medium  and 
their  validity  is  not  limited  to  elastic  shells  alone. 

Finally,  we  briefly  describe  here  the  results  of  a recent  study  by  Naghdi 
and  Tang  f lB  J concerning  control  Inbl  o deformations  that  can  bo  maintained,  in 
the  absence  of  body  force,  in  every  isotropic  elastic  membrane  by  the  applica- 
tion of  edge  loads  and/or  uniform  normal  surface  loads  on  the  major  surfaces  of 
the  thin  shell-like  body.  The  static  solutions  of  finitely  deformed  membranes, 
which  are  valid  for  both  compressible  and  incompressible  materials,  are  obtained 
with  the  use  of  a strain  energy  response  function  which  depends  on  the  metric 
tensor  of  the  membrane  in  its  deformed  configuration.  The  main  results  are 
summarized  by  several  theorems  and  their  corollaries  in  accordance  with  throe 
mutually  exclusive  cases  for  which  the  initial  undeformed  surface  of  the  membrane 
(which  may  be  a sector  of  a complete  or  closed  surface)  is,  respectively,  develop- 
able, spherical  and  a surface  of  variable  Gaussian  curvature  satisfying  certain 
differential  criteria.  The  corresponding  deformed  surfaces  are,  respectively, 
a plan  or  a right  circular  cylinder,  a sphere  and  a surface  of  constant  mean 
curvature.  These  results  are  exhaustive  in  that  they  represent  all  finite 
<b  formation  solutions  possible  in  every  isotropic  elastic  material  characterized 
oy  the  strain  energy  response  mentioned  above.  Also  discussed  in  [l8J  are  same 
special  cases  of  the  general  results  and  several  families  of  solutions  in  terms 
of  an  alternative  description  which  should  be  useful  in  application  and  which 
permit  easy  interpretations. 
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. . . L.  \n  t • menl  of  1 shell -lik<  bo  ly  with  uni  form  thickness  it:  reference 
configuration  showing  the  middle  surface  0 and  th*  major  surface; 

* = i h/2.  Also  indicated  art  the  position  vector  H from  n fixed  origin 
to  a point  on  the  middle  surface  £ = 0,  the  position  vector  i to  a point 
of  the  region  of  space  occupied  by  the  shell  in  the  reference  configura- 
tion and  the  unit  normal  A0  to  the  middle  surface. 
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